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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE 23^ 

THE EFFECT OF AN ARBITRARY SURFACE-TEMPERATURE VARIATION 
ALONG A FLAT PLATE ON THE CONVECTIVE HEAT TRANSFER 
IN AN INCOMPRESSIBLE TURBULENT BOUNDARY LAYER 
By Morris W. Rube sin 


SUMMARY' 


An analysis is performed to determine the effects on the rates of 
convective heat transfer produced by variations of the surface tempera- 
ture with distance along a flat plate. The analysis is confined to the 
case of a low— speed turbulent boundary layer on a flat plate oriented 
parallel to the free— stream velocity. All the physical properties of 
the fluid are considered to be constant. The temperature distribution 
in the boundary layer and the local heat— transfer rate at the surface 
are obtained in the form of definite integrals, with the integrands con- 
taining the prescribed surface— temperature variation or its derivative 
with respect to the distance along the plate. Numerical evaluation of 
these integrals permits application of the results to any type of sur- 
face temperature distribution that may occur physically. 

The basic solution of this analysis, one applying to a stepwise- 
discontinuous surface temperature in which the temperature is constant 
on either side of the discontinuity, is correlated with a set of exist- 
ing experimental data. This work is extended to show that the existence 
of a st epwise-discont inaigus^aur^&e^temP.eratm'aiUP.s.t rea m of a plug-type 

hea t meter makes the resu lts,,,. obtained, from thi s type of an instrument 

subject to large .errors . 

For convenience, algebraic equations are presented for the surface 
temperature , the surface temperature gradient with respect to the dis- 
tance along the plate, and the local convective heat— transfer rate. 

These equations are in the form of power series, with positive exponents, 
in terms of the distance along the plate. When one of the quantities, 
surface temperature or local heat transfer, is specified and the coef- 
ficients and exponents of the particular power series are evaluated, 
the other quantities, together with the total heat— transfer rate, can 
be determined directly by substituting these coefficients and exponents 
into the respective power series. 
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INTRODUCTION 


Several experiments and analyses have teen performed which indicate 
that a variable surface— temperature distribution can produce either a 
marked increase or decrease in the local and average convective heat- 
transfer rates to or from a surface. In reference 1 an approximate anal- 
ysis is described in which the effects on the convective heat— transfer 
rates produced by a stepwise-discontinuous surface temperature on a flat 
plate were determined. These effects were determined from solutions of 
the integral form of the boundary— layer-energy equations for both lami- 
nar and turbulent boundary layers. Further, these solutions were 
restricted to the assumption of constant physical properties for the 
fluids. It was found that extremely large effects occur in the region 
directly downstream of the discontinuity in the surface temperature. 

In reference 2, there is described an analysis in which the analyt- 
ical solution of the boundary— layer equations was determined for the 
compressible, laminar boundary layer on a flat plate. The surface tem- 
perature in this analysis was represented by an arbitrary polynomial in 
terms of distance along the plate. It was found that even a continuous 
surface— temperature variation could produce large effects on the convec- 
tive heat transfer in the laminar boundary layer. 

In reference 3 > there is described an analysis in which the solu- 
tion of the boundary— layer energy equation was determined for laminar 
flow on the surface of a wedge. When the velocity distribution in the 
boundary layer was assumed to be linear with the distance normal to the 
surface and all the physical properties of the fluid were assumed to be 
constant, a solution of the energy equation was. obtained analytically 
for a stepwise-discontinuous surface temperature. For the case of a 
flat plate this solution gives results in agreement with the results of 
the solution for the laminar boundary layer in reference 1. This solu- 
tion (reference 3) for the case of a stepwise-discontinuous surface tem- 
perature was then extended to an arbitrary surface temperature through 
the employment of an integral solution similar in idea to that of 
Duhamel (reference 4). Recently another analysis (reference 5) was per- 
formed determining the heat transfer from a body on which the velocity 
at the edge of the boundary layer aid the surface temperature vary. 

This analysis required the same limiting assumptions as in reference 3 
although the mathematical details were somewhat different. Again it was 
found from these analyses that the surface— temperature variation influ- 
enced the convective heat— transfer rates considerably. 

It has also been shown experimentally that a variation in the sur- 
face temperature with distance along the surface produces large effects 
on the convective heat transfer. In reference 6, an experiment is 
described in which the effect of unheated starting sections on the aver- 
age heat transfer was determined from a cylindrical probe having a 
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constant surface temperature in the heated region and with, its axis par- 
allel to the free— stream velocity. It was found that the length of the- 
uhheated starting .section does influence the average rate of heat trans- 
fer on the remainder of the probe. 

In the experiment described in reference 1 , it was found that val- 
ues of the local heat— transfer rates measured on a cone where severe 
surface— temperature variations occurred deviated from predictions based 
on theories in which the surface temperature is assumed constant. These 
results were obtained for both laminar and turbulent boundary layers. 

Recently an experiment was performed to determine the effect of a 
stepwise-discontinuous surface temperature on the local heat transfer 
in the turbulent boundary layer of a flat plate, reference 8. 1 These 
data show a marked effect on the local heat transfer, and will be cor- 
related with the results of the analysis of this paper. 

A quantitative determination, either experimental or analytical, 
is necessary to reveal the effect of a continuous variation of surface 
temperature with distance on the convective heat transfer in the turbu- 
lent boundary layer. It is the purpose of this paper, therefore, to 
make this determination in an approximate analytical manner subject to 
the limiting assumptions of a flat plate, of constant physical proper- 
ties of the fluid, and of no frictional dissipation of energy within 
the boundary layer. 


STMBOLS 

coefficients of power series, °F per foot -11 

39+a8m 

symbol defined by equation (22), foot 140 


B 
C 
c 

P 
F 

i~ 

The author wishes to acknowledge his indebtedness to Mr. Steve Scesa 
and the Department of Mechanical Engineering, University of California, 
Berkeley, for their kind permission to use hitherto unpublished data 
in this report. 


coefficients of power series, dimensionless 

. . , —1/35 

symbol defined by equation (42), foot 
coefficient defined by equation (25), dimensionless 
specific heat at constant pressure, Btu per slug, °F 
coefficient in equation (37 ) y dimensionless 
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G coefficient in equation (37 ) > dimensionless 

h local heat-transfer coefficient, Btu per second, square foot, °F 

h* local heat— transfer coefficient defined hy equation (66), Btu. per 

second, square foot, °F 

h average heat— transfer coefficient, Btu per second, square foot, °F 
H coefficient defined hy equation (40), dimensionless 
k summation index, dimensionless 

1 distance along flat plate to point of the stepwise discontinuity of 
the surface temperature, feet 

m parameter introduced in equation (20), dimensionless 

n summation index, dimensionless 

Nu Nusselt number (hx/k), dimensionless 

P function defined hy equation (34), dimensionless 

Pr Prandtl number (pcp/x), dimensionless 

q local heat— transfer rate per unit area, Btu per second, square foot 

Q total heat— transfer rate per unit width, Btu per second, foot 

Be Reynolds number (u Q x/v), dimensionless 
s total length of heated portion of plate, feet 
t temperature of fluid in the boundary layer, °F 

u local velocity in the boundary layer parallel to plate, feet per 

second 

u q free-stream velocity, feet per second 

U upper limit of integral defined in equation (48), feet 

v local velocity in the boundary layer normal to plate, feet per 

second 

W distance from leading edge of flat plate to rear of heated section 
following a surface— temperature discontinuity, feet 
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x distance along plate from leading edge, feet 
y distance normal to plate surface, feet 
Y n function defined "by equation ( 65 ) 

z temperature ratio defined by equation ( 38 ), dimensionless 
5 thickness of the flow boundary layer, feet 
A thickness of the thermal boundary layer, feet 

9 exact temperature distribution in the boundary layer for a plate 

at constant temperature preceded by an unheated starting section, 
dimensionless 

9 * approximate temperature distribution in the boundary layer for a 
plate at constant temperature preceded by an unheated starting 
section, dimensionless 

X thermal conductivity of fluid, Btu per second, square foot, °F per 
foot 

p absolute viscosity of fluid, pound-seconds per square foot 
V kinematic viscosity of fluid, square feet per second . 

2 4 

p mass density of fluid, pound-HBeeonds per foot 
t local shear stress, pounds per square foot 

Subscripts 

e effective property for turbulent flow 
o free— stream condition 

w surface condition 

p referring to surface of heat-meter plug 

s referring to surface of surrounding material 



6 


ITACA TN 23^5 


AHALISIS 

General Expressions for the Local Heat— Transfer Rate and 
the Temperature Distribution in the Boundary Layer 


The method of analysis in this paper will be to determine an approx- 
imate solution of the boundary— layer energy equation in integral form for 
a stepwise-discontinuous surface temperature , and then to extend this 
solution to apply to an arbitrary surface temperature by employing an 
integral solution similar in idea to that of Duhamel (reference 4) . 

Boundary— layer temperature distribution.— The energy equation for 
the turbulent boundary layer on a flat plate, neglecting viscous dissi- 
pation, can be expressed as 


pc^ f\L ^ + v "\ 

P V a 7 ) 



(1) 


where X e , p, and Cp are assumed to be independent of temperature. 
Suppose an exact solution of this equation exists for the case of a 
plate which is unheated and at the free— stream temperature t Q in the 
region x < L, and is heated and at the temperature t w in the region 
x > L. For the region x > L 


12s - * e>(x,y,L) 

■to-^w 

where 

0(x, 0,L) = 0 
0(Lf,y,L) = 1 
e(x,°o,L) = 1 


(2) 


Because of the linearity in temperature t of equation (l), various 
solutions may be added to satisfy desired boundary conditions. For 
instance, the temperature in the boundary layer on a flat plate which 
is maintained at t Wi , different from t Q , for x < L and t^ , 
different from t Q , for x > L is given by 

t = tw 2 +(t o -t Wl )0(x,y,O)+(t Wi -t W2 )0(x,y,L) 


(3) 
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in the region x > L. Each term' of equation ( 3 ) satisfies the energy- 
equation and the sum of the terms satisfies the desired "boundary condi- 
tions ; thus, the temperature represented by equation (3) can be con- 
sidered as the formal solution of the problem. The idea underlying equa- 
tion (3) may be extended to a large number of temperature discontinuities 
^7 



k=i 




(*) 


where x > L, . Equation (4) may be rewritten as 

K — X 


n 


t = t o -.0(x,y,O) + ^ t„ k - eCx^Lk-J 

k=i 


(5) 


Equation (5) is restricted to x = L^. If n is allowed to increase 
indefinitely, the interval (L^— ) becomes small and equation (5) 

may be rewritten in integral form as 


t(x,y) 


t Q 0(x,y,o) 


■+ 



t(L,a) — dL 


(6) 


Equation (6) constitutes the general expression for the temperature dis- 
tribution in the boundary layer for an arbitrary surface— temperature dis- 
tribution. 


Local heat— transfer rate.— The convective heat— transfer rate per 
unit area from a surface, where X = X e , is given by 


q = — X 


St 



(7) 


In order to introduce a temperature potential in the heat-transfer expres- 
sion, equation (6) is integrated by parts which yields 

^x 

t(x,y) ■* t(x,0) + [t o -t(0,0)] 9(x,y,0) - / ^M Vdl 

J Sl 

O 


(8) 
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after use of the boundary conditions of equation (2). The local rate 
of heat transfer obtained by differentiating and substituting equa- 
tion (8) into equation (7) is 

q(x) = -x[t Q — t(o,o)] deU,y,°) | + \ r fo .. ( . L /0) . be I ^ ( 9 ) 

* y=o O y=o 


When the local heat— transfer coefficient is defined by 

h(x,L) = X, d e ( x j7,L) (10) 

y=0 

equation (9) becomes 

x 

q(x) = h(x,0) [t(0,0)-t o 3 + / h(x,L) dL (ll) 

J 8l 

o 

Equation (11) constitutes the general expression for the local heat- 
transfer rate per unit area. Note that h(x,0) represents the heat- 
transfer coefficient on a plate having a constant surface temperature. 

It is interesting to note that the local 'eat— transfer rate per unit 
area is comprised of the local heat— transfer rate on a plate of constant 
temperature, equal to the leading— edge temperature of the plate, plus 
corrective terms introduced by the variation of the surface temperature. . 


Approximate Solution for a Stepwise- 
Discontinuous Surface Temperature 


General solution. — The purpose of this section is to obtain the 
value of h(x,L) to be used in equation (ll). Because of the complex- 
ity of the problem, an exact solution 0, or indirectly h(x,L), cannot 
be determined at present and an approximate approach is indicated. Thus, 
the a na lysis presented in reference 1 will be repeated herein in a more 
general fashion. 

The equation which will be solved is the boundary— layer energy 
equation A 

v _ pc P d r , 

“(t-t 0 )dy (12) 

o 
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The method of solving this equation is to transform it. to an ordinary 
differential equation in terms of x as the independent variable 
and A as the dependent variable. The resulting differential equation 
is then integrated. 

Suppose the velocity term u in the right member of equation (12) 
is given by 


u _ 


u o 



(13) 


where 


5 « 0.37 u 


- 1/5 1/5 4/5 

O V X . 


(14) 


These are the well-known l/7— power— law relationships (reference 9). In 
addition, suppose 


6 * 



(15) 


even though 5 and A may not be of the same order of magnitude . When 
the velocity and temperature expressed by equations (13) and (15) are 
substituted into the right member of equation (12), there is obtained 
after simplification 


h 

X 



(16) 


when A is less than 5. It should be noted that for the case of a 
constant plate temperature in air A equals 1.26. When the integra- 
tion of equation (12) is performed with these limits on the thicknesses 
of the boundary layers it is found that the results of the heat— transfer 
coefficient differ from those of equation (l 6) by the order of 0.1 per- 
cent. Therefore, even though A is greater than 6 for some cases, 
the results of the following analysis are believed to be valid over tlie 
entire plate to the order of 0.1 percent for air. 

In order to transform equation (16) to an ordinary differential . 
equation, it is necessary to express the left member of equation (l6) in 
terms of the local characteristics of the boundary layer. Many alterna- 
tive expressions may be obtained in this manner. This avoids imposing 
a prescribed variation of the heat— transfer coefficient h in terms 
of the distance along the plate x, which would make the problem 
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trivial. To determine which of the alternative solutions is applicable, 
correlation of a general solution depending on a single parameter will 
be made with the data of reference 8 . The value of the parameter which 
achieves the best correlation of the data and the analysis will then be 
used in the remainder of the report. 

In order to obtain these alternative expressions for h in terms 
of local characteristics, it is necessary to restrict the analysis, for 
the moment, to the case of a constant surface temperature. For the case 
of a constant surface temperature, the Colburn analogy between local 
heat— transfer rate and skin friction (reference 10 ) is 


^-(Pr ) 2/3 --4 (17) 

P“o°p P u o 

The right member of equation ( 17 ) is expressed, from 1 / 7 — power-law con- 
siderations, as 


= 0.0225 u“ 1/4 V 1/4 B“ 1/4 (18) 

pu Q 2 

Further, it is shown in appendix A that for the case of a constant sur- 
face temperature the ratio of the thickness of the thermal boundary 
layer to the flow boundary layer is given by 

A — 7/12 

| - (Fr) (19) 

o 

When equations ( 17 ), (l 8 ), and (19) are combined so as to express h/X 
in terms of Pr, 5 , and A, the variables Which appear in equation ( 16 ), 
there results 

/ n 3/4 (s— 29m) — i+4m 

^ = 0.0225 Pr 48 A 4 5 m (20) 

The thermal boundary— layer thickness was introduced into this equation 
by raising both sides of equation ( 19 ) to the m^-th power. Equation (20) 
constitutes a general expression for the local heat— transfer coefficient 
in terms of the local boundary— layer characteristics and the Prandtl 
number. Although this equation was derived for heat transfer from a 
plate at a constant surface temperature, it will be assumed to apply to 
the case of a flat plate with the region x < L unheated and at the 
free-stream temperature. It is interesting to note that there are 
limits on the possible variation of the parameter m. For instance. 
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in equation (20) when m is less than — l/4 it is found that the heat- 
transfer coefficient increases with, an increasing thermal boundary— layer 
thickness. When the parameter m is greater than 9/28, it is found 
that the heat— transfer coefficient increases with a decreasing Prandtl 
number. These variations are in contradiction with existing knowledge 
concerning convective heat transfer. The values of the parameter m 
must , therefore , lie in the interval —l/k < m < 9/28. 

When equations (16) and (20) are equated, there results after sim- 
plification 


dA 

dx 


— A = Ax 

lOx 


gam+4 

33 


A 


n-rggm 

28 


( 21 ) 


where 


7 m+i y \ _ 39H-g8m _ 39+28m 

A = (0.2025)(0.37) ■ 7 0^ ) 140 £r 48 (22) 

Equation (21) is of the Bernoulli type which can be integrated to yield 


28 


r 

/ 224D1+312 19031+273 39+2031 v 


39+28S1 


Jr A 

{ x 280 — L 280 x 280 

)] 


(23) 


when the boundary condition A = 0 at x = L is imposed. Equation (23), 
therefore, gives the thermal boundary— layer thickness on a plate which is 
unheated and at the free-stream temperature t Q for x < L and is 
heated and at the temperature t w for x > L. When equation (23) is 
substituted into equation (20) there results 


h(x,L)x 

\ 


= C(m) 






28HH-39 

40 


7+28m 

39+88m 


(2k) 


The coefficient in equation (24) is given by 

3g gam — 1 7+28m 

. , . . , 39 + 28 m . . 33+28m . „ ,39+2 8m . 

C(m) = (0.0225) (0.37) (0.0778) (25) 


In the range — l/4 < m < 9/28 the value of C(m) is essentially con- 
stant and equal to 0.0288. It is apparent from equation (24) that when 
m is greater than — l/4 and when x approaches L the effect is to 
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increase the local heat— transfer coefficient many fold. When L equals 
zero this equation degenerates to the customary heat— transfer expression 
based on the 1/7— power laws for velocity and temperature distribution 
(reference 11). 

Correlation with experiment.— To determine which value of m makes 
equation (24) more nearly conform to the measurements, comparison is made 
of this solution with the data of Scesa in reference 8. The data of 
reference 8 were determined in a well— eontr oiled experiment on a flat 
plate oriented parallel to the free-stream velocity. Free-stream speeds 
of 50 to 70 feet per second were employed. The plate had an unheated 
starting section which was followed by a heated section in which there 
was an additional temperature discontinuity. All the basic data except 
those from heater elements which can be considered guard heaters are 
plotted in figure 1, It is observed that the data lie consistently above 
the theoretical line for a plate at a constant temperature. 

The temperatures in the heated region of the plate are defined as 
t Wl and upstream and downstream from the temperature discontinuity, 

and the temperature of the unheated starting section is t Q , the free- 
stream temperature . The length of the unheated starting section is Li, 
and the distance from the leading edge of the plate to the temperature 
discontinuity in the heated portion of the plate is L 2 . The local heat- 
transfer rates per unit area are then given, respectively, by 


q = h(x,L 1 )(tw 1 — to) 


(26) 


for Li < x < L 2 , and 


q = h(x,L 1 )(t Wi -t 0 ) + h(x,L 2 )(t W2 -t Wl ) (27) 


for x > L 2 . The heat— transfer coefficients in this experiment were 
defined in terms of the local temperature difference between the surface 
and the free stream. The Nusselt numbers corresponding to those meas- 
ured in each of these regions are given by equation (24) as 


L=0 


28DH-39 



7+28 m 
39+2 8m 


Hu = (Nu). 


(28) 



NACA TN 23^5 


13 


and 


Nu = (Nu) 

L=0 


~^w i~~k o 
tw 2 — to 




28m+39_— 
40 


7+28m 
39+28 m 


+ 


twp-t- 


W.l 


t w 2 -to 



> v 28Ittt-39 

7+2 8m | 

- 

39+28 Ht l 

l-( 

hg j 4 0 



/ 

J 


(29) 


If the general solution, equation (24), yields results which conform to 
experiment, then the experimental Nusselt numbers divided hy the brack- 
eted quantities should yield experimental Nusselt numbers which corre- 
spond to (Nu)^_q. These adjusted data are shown for two values of the 

parameter m, m = 0, and m = 9/28, in figures 2 and 3> respectively. 
It should be noted that, m = — 1/4 corresponds to the data of figure 1, 
that is , the bracketed quantities are equal to unity. The value m = 0 
apparently correlates the data very well. This corresponds to the 
results of reference 1, that is, the local heat— transfer coefficient is 
dependent on the thermal boundary— layer thickness and the Prandtl number 
and is independent of the flow boundary— layer thickness. For m = 0 
equation (24) becomes 


h(x,L)x _ 0 < o 288 
X 



(30) 


Equation (30) will be used in the remainder of the analysis. 

Average heat-transfer coefficient determined from particular 
solution.— In many cases it is desirable to know the average heat- 
transfer coefficient over a region at constant temperature preceded by 
an unheated starting section. The total heat transferred in the region 
per unit width is given by 




h(x,L)(t w -t 0 ) 


&x 


(31) 


The symbol ¥ represents the distance from the leading edge of the 
plate to the rear of the region considered. The average heat-transfer 
coefficient over this region is defined as 


h = 


Q 


(tw-t 0 )(¥-L) 


(32) 
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When equations (30), (31), and- (32) are combined there results 


4/5 


h(W,L) = h(W,0 ) 


40 (l/w) 

39 i-(l/w) 


(w) 


(33) 


where 


P 




T,- 71/3S (l-,)- 7 / 33 jfa, 


(3<t) 


Equation (34) was evaluated numerically and the results are shown in 
figure 4. 

The experimental data obtained by Jakob and Dow (reference 6) on 
the average heat transferred from a probe in axial air flow preceded by 
unheated starting sections corresponds, in form, to the results obtained 
from equations (33) and (34). It cannot be assumed tacitly, however, 
that the data obtained from probes should compare identically with 
theory based on a flat plate. In particular, the data obtained on the 
probe for a constant surface temperature were about 13 percent lower 
than corresponding theoretical or experimental results on a flat plate. 
As any theoretical comparison of probes and flat plates indicates that 
the rate of heat transfer from the probe should be the greater, there 
appears to be some unexplained reason for the low values of the data. 

In view of this, no extensive comparison of the data and the results of 
the present analysis are made. A s ingle_c ompar i s on for the case of an 
unheated starting section reveals that h/h(W,0) = 1.37 from the Jakob 
and Dow results , while the corresponding result from the present analy- 
sis is h/h(W,0) = 1.45. This comparison is favorable when considera- 
tion is made of the questions expressed previously concerning the com- 
parison of results from probes and flat plates. 


Particular solution applied to general single surface— temperature 
discontinuity.— When the region of the -plate preceding the temperature 
discontinuity is at a temperature t Wi other than the free-stream temr- 

perature, the local and average heat— transfer coefficients in the region 
L < x < W at a temperature are given by 


1 t w J -t o 

’ fc w 2 — fc o ^w 2 ~^o 




) 


39/40 f 7 / 39 


h(x,L) = h(x,0) 


(35) 
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h(W,L) = h(W, 0 ) 


5 Li-(l/w) 

^ i-(l/w) 


4/5 


^Wx^o + to (L/W ) 4/5 P(L/W) V^x I 
tw 2 -to + 39 i-(l/w) tw 2 -t 0 


( 36 ) 


For certain applications in a small region (W— L) it is interesting to 
know the ratio of the average heat— transfer coefficient to the local 
heat— transfer coefficient which would prevail at the center of the region 
if the plate were at a constant temperature. This ratio is obtained by 



where F(l/w ) and G(l/w) were evaluated numerically. If 


z 


^*w 2 % x 
t O 


( 38 ) 


equation ( 37 ) may be rewritten as 


where 


h(W,L) 



= F(L /¥) + H(L/W)z 


h(l/w) « G(l/W) - f(l/w) 


( 39 ) 


(to) 


The numerical values of F(l/W) and l(L/¥) are plotted in figure 5 . 


Extension of Approximate Solution to the Problem 
of an Arbitrary Surface Temperature 


Boundary— layer temperature distribution. — The approximate temper- 
ature distribution determined from equations (15) and (23) with m = 0 
is 
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e* = 


0 


By j — 

3X2/280 273/280 39/280 \41/39 


- L 


X 


/ 


where 


B = 1.15 u D 1/35 v 1/35 Pr 1/12 


(to) 


(to) 


At the outer edge of the boundary layer 0* becomes equal to unity; 
therefore, equation (to) is valid only in the region 


312/280 273/280 39/280\28/39 

y< , £J -L ». ) 



(43) 


For values of y greater than given by the inequality (43) 

0 * = 1 

For the range of y indicated by the inequality (43) 

. _l/7 39/280 T — 1/40 

d0» _ B y * -L 

<3L 10 f 312/280 273/280 39/280\ 43/39 

lx — L x J 

For greater y 


Ml 


= 0 


(to) 


(to) 


(46) 


When equations (4l), (45), and (46) are substituted into equation (6) 
there results 


t(x,y) = By 


1/7 


v . + .i r u . 

, , 39/280 —1/40 

t(L,0) x L ' dL 

4/35 10 J 

L z o 

/ 312/280 __ L 273/280 ^.39/280^ 43/39 


( 47 ) 
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where 


g _ ^.273/280 _ b 39/4 y 39/28 x -39/280^ 


x 280/273 

') 


IT 

(48) 


1/7 4/35 .. 

When By = x it is apparent from equation (4o) that U = 0, and 

consequently from equation (4-7) t (x,y) = t Q . It is not obvious, however, 
that the left member of equation (4-7) equals the prescribed surface tem- 
perature at y = 0, It is shown in appendix B that equation (47) does 
satisfy the boundary conditions at y = 0. 

Local heat— transfer date.— The general expression for the local 
heat— transfer rate is obtained by substituting the local heat— transfer 
coefficient given by equation (30) into equation (ll) 


q = h(x,0) 


r /■_ \ 39/4o n 

■4 

—7/39 

H 

1 

dL 

<■ 


(49) 


For an arbitrary surface— temperature distribution it is necessary 
to evaluate the integrals in equations (47) and (49) numerically. There 
is, however, a type of surface distribution for which the local heat- 
transfer rate can be determined easily by analytical treatment. This 
will be shown in the next section. 

Surface— temperature distribution represented by -power series,— Let 
the surface temperature be represented by 


(50) 


t(x,0) - t Q = 2^ a^x 11 n > 0 
n 


where n is not necessarily an integer. When equation (50) is differ- 
entiated with respect to x there results 


= y 

8x 


(51) 


The substitution of equation (51) at x = L into equation (49) yields 


q = 


h(x,0) ‘ a D + J ^ na^L* 1 * 

{ o n 


l-l * 

.Xy 


39/40 


-7/39 


dL 


(52) 


This equation may be rewritten as 


q - h(x,0) 


a o + 7 m n 


-** [* -© 


39/40 —7/39 


dL 


(53) 
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. , .39/40 

Letting a = (L/xj allows transformation of the integral terms of 

equation (53) so that the equation reads 


q = h(x,0) 





±9n- 


a’ 


39 


n o 


— T/39 

(I-a) da 


(5*0 


The integral terms of this equation are Eulerian integrals of the first 
kind (reference 12), the values of which are given by the Beta function 
which may be expressed in Gamma functions as 



where n is greater than or equal to zero. The values of the Gamma 
functions may be obtained from tables in reference 13 . In view of equa- 
tion (55 ) f equation (5*0 becomes 



When the recursion equation 

r(y+i) = yr(y) (57) 


is used, equation ( 5 6 ) may be rewritten as 


1 = 



n 



(58) 


It is apparent that once the type of surface— temperature distribution 
is specified, that is, values of a n specified, the distribution of 
the local heat— transfer rate per unit area can be determined directly 
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from equation (58). The total rate of heat transfer per unit width Q 
is obtained by integrating equation (58) from x = 0 to x = s 


Q = h(s,0) 




n ^ 


S n+1 



(59) 


The term h (s,0) is the average heat— transfer coefficient for the case 
of a constant surface temperature determined from the expression 


= 0.036^^-^ Pr 1/3 


(60) 


For convenience, equations (5I) , (52), (58), and (59) are put in 
the dimensionless form 


t(x,0) 


= 1 + 


I <*)' 


t(x 


£) N 

to / 


n— 1 


fir- ?■ 


7, • I *•'■(;)* 


t 0 h(x 


(61) 


(62) 


(63) 


q = y ^ t 

t o 'sK(s,0) L d n+1 _) , » 


m 


In these equations 


Y = 


ho 

r ( ~ n+l 
39 


Mi) 


r[ + 


39 39 


32 \ 
qq / 


(65) 
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The term. Y n is plotted as a function of n in figure 6. If the heat- 
transfer coefficient for the case of a variable surface temperature is 
defined in terms of the local temperature difference 


* 


h 


q 

t(x,0) — t Q 


(66) 


then 



(67) 


Inspection of the right member of equation (67) reveals several inter- 
esting facts. For a problem in which the denominator has a local value 
of zero, the numerator will be finite in all but a very exceptional case; 
thus h*/h(x,0) will be infinite at this point. For another case the 
numerator may be zero locally while the denominator will probably be 
finite; thus h*/h(x,0) will be zero at this point. It is seen, there- 
fore, that h*/h(x,0) can attain any numerical value, depending on the 
problem considered. 


In summary, it is noted from equations (6l) , (62), and (63) that the 
values of the coefficient A n may be determined by a prescribed distri- 
bution of surface temperature, a prescribed distribution of the surface 
temperature gradient, or a prescribed distribution of the local heat- 
transfer rate per unit area. Once the coefficients A n are evaluated 
from the prescribed quantity the other two quantities, together with the 
total heat— transfer rate, can be determined directly. 


EXAMPLES 

Plug— Type Heat Meter 


Convective heat— transfer rates at a surface are often measured con- 
veniently by determining the heat— transfer rate to or from a plug which 
is thermally isolated from the surrounding material. The plug can be 
used as a steady-state device in which heat is generated electrically 
and the temperature of the plug is maintained constant. The usual prac- 
tice, however, is to use the plug as a transient device by measuring, the 
temperature— time history of the plug. The local heat— trans f er rate 
through the surface of the plug is represented by the product of the 
rate of change of temperature and the thermal capacity of the plug. 
Kegardless of the convenience of using such a device, the present 
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analysis indicates that the measurements are subject to certain inherent 
errors in addition to those which may arise from imperfect thermal iso- 
lation. 

It is not possible to design a simple plug— type heat meter which 
will always maintain a temperature identical to the temperature of the 
surrounding material. If an extremely thin insulating material is used 
between the plug and the surrounding material, a surface— temperature dis- 
continuity will occur at the seam of the plug and the surrounding mate- 
rial. If the instantaneous surface temperature of the plug is defined 
as tp, while the surface temperature of the surrounding material to the 
leading edge of the plate is denoted t g , it is apparent that equations 
( 36 ) through (4o) are applicable to this problem, where tp = ty 2 and 
= twx- 

It is desirable to know how well a plug— type heat meter measures 
the local heat— transfer coefficient which would exist at the position of 
the center of the plug if the plate were at a constant temperature. 
Equation (37) is directly applicable to this problem. For example, sup- 
pose the instantaneous temperatures of a heat-meter plug, the surrounding 
material, and the free stream are 125° F, 100° F, and 200° F, respec- 
tively. The dimensions of the plug are L = 2k inches and W = 25 inches. 
It is found from equation ( 38 ) that z = -0.33> and- I*/w = O. 96 . From 
equation ( 39 ) and figure 5 


- h( 25 , 2 lQ 

h(24.5,o) 


( 68 ) 


Thus, the plug— type heat meter measures a heat— transfer coefficient which 
is only 62 percent of the local heat— transfer coefficient which would 
have existed at the position of the center of the plug had the plate been 
at a constant surface temperature. 


Although this example is arbitrary, the conditions are by no means 
implausible. It should be noted from equation (39) that more severe 
surface— temperature conditions can result in negative values of, . 


h(W,L) 
h (¥+ 1 / 2 , 0 )' 


which means that measurement of' heat transfer by a plug 


would be in a direction opposite to that which would normally prevail on 
a plate having a constant surface temperature. These results indicate 
that extreme caution should be exercised in the interpretation of the 
data obtained by plug— type heat meters. 
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Constant Heat— Trans f er Hate 


Suppose it is desired to maintain a surface— temperature distribu- 
tion. such that the local heat— transfer rate per unit area is constant 
along a plate. From equation (30) it can be seen that h(x,0) is 
inversely proportional to the l/5— power of x. Equation (63) can there- 
fore be written as 


q = t D h(s,0) 1/5 £ ) <«9) 


It is apparent that to maintain q constant, n must have the value 
l/5 only. Equation (69) becomes 


q = t Q h(s,0) A 1/5 Y 1/5 


(TO) 


Therefore, 


1/5 t Q h(s,0) Y 1/5 


(71) 


The required surface— temperature variation as given by equation (6l) is 


t(x,0) 


1 + 


q 

t Q h(s,0) Y 1/5 



(72) 


Further 


* 

h 


h(x,0) 



1.06 


(73) 


This latter equation indicates, from the values of figure 6, that the 
local heat— transfer coefficient on a plate having a constant heat- 
transfer rate differs from that on a plate at a constant temperature 
by only 6 percent. 
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Surface Temperature Represented "by Polynomial 
Suppose the temperature distribution on a flat plate is given hy 

* = I* 5 + °- 5 (j) W 

Then in equation (.61)," n = 0 and n = 1, and A Q = 0.5 and Ai = 0.5. 
From equation ( 63 ) and figure 6, the local heat-transfer rates are given 

"by 


2— = 0.5 + 0.612 

t 0 h(x,0) 



(75) 


The ratio of the local heat— transfer coefficient to the heat— transfer 
coefficient which would exist if the plate were at a constant surface 
temperature is 


h* _ 0.5 + 0 . 612 (x/s) (6) 

h(x,0) 0.5 + 0.5 (x/s) 

The values of the members of equations (7*0 and (75) are plotted in 
figure 7 . It can be observed that although the temperature -variation 
along the plate increased the over-all temperature difference between 
the plate and the free stream by 100 percent, the local heat— transfer 
coefficient differs from that on a constant— temperature plate by a maxi- 
mum of about 11 percent. 


CONCLUSIONS 


The conclusions of this report are subject to certain limitations 
inherent in the assumptions of the analysis. The following assumptions 
have been used: 

1. All physical properties of the fluid are constant. 

2. Frictional dissipation of energy within the boundary layer is 
negligible. 

3. The velocity distribution in the boundary layer is of the 
l/7— power form. 
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The temperature distribution in the boundary layer is of the 
l/7— power form even when the thermal boundary— layer thickness differs 
considerably from the flow boundary— layer thickness. 

5. The local heat— transfer coefficient determined on a plate having 
a constant surface temperature applies to a plate having a variable sur- 
face temperature when it is expressed by an equation based on the local 
flow and thermal boundary-layer thicknesses instead of the distance along 
the plate. 

Although the analysis based on these assumptions is correlated with 
experimental data, this is not a verification of the individual assump- 
tions of the analysis, but rather, is a justification of the use of the 
end results. 

The foregoing analysis, together with the examples cited, has indi- 
cated that a variation in the surface temperature with distance along a 
flat plate influences the local convective heat transfer to an extent 
which depends on the type of variation. In general, a sudden surface- 
temperature jump, or discontinuity, produces extremely large increases 
or decreases in the convective heat transfer directly downstream of the 
position of the discontinuity. From this it can be concluded that meas- 
urements by a heat meter of the plug type can, in most instances, deviate 
considerably from the heat transfer which would normally exist at the 
location, of the instrument. A continuously variable surface temperature 
was shown, in general, to have a smaller effect on the convective heat 
transfer than a sudden temperature discontinuity. For instance, it was 
shown that the local heat— transfer coefficient on a plate having a con- 
stant heat— transfer rate along the surface, with the necessary l/5— power 
variation of the surface temperature, differs from that on a plate with 
a constant surface temperature by only 6 percent. For the case of con- 
tinuously varying surface temperatures , large effects on the local heat- 
transfer coefficient are expected only where the surface temperature 
approaches that of the free stream at a point other than the leading 
edge of the plate. 


Ames Aeronautical laboratory. 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Feb. 9, 1951, 
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APPENDIX A 

RATIO OF BOUNDAET-LAYEE THICKNESSES ON PLATE 
WITH CONSTANT SURFACE TEMPERATURE 


If the 1 / 7 — power— law velocity distribution 


it ■ (!)‘ 


is substituted into the von Karman momentum integral for a flat plate 


— / u(u 0 ^i)dy 

P dx J 
0 


there results 


t _ d 
P dx 


<sn>->fcr> 


When the indicated operations are performed, there results 


7 dS 


pu 0 2 72 dx 


Similarly, if the 1 / 7 - power temperature distribution 


ttf-ty \ A 


and equation (Al) are inserted into the integral form of the energy 
equation 


h = pCp — f u(t-t Q ) dy 
, V^o dx J 


there results 


h d ( 7 A 8/7 ' 

P ° dX\ 72 gl/7 , 
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when t w is constant and it is postulated that 5 is greater than A. 

For a constant surface temperature the Colburn modification of the 
Reynolds analogy is expressed as 


pc p u Q 


(Pr) 


2/3 


m 0 s 


(A8) 


When equations (kk) and (AT) are substituted into equation (A8) , there 
is obtained 


_d_ 

dx 


2/3 

(Pr) 


a 6/7 1 
# /7 



(A9) 


Therefore 


2/3 

(Pr) 


a 8/7 

6 1 / 7 


8 + C 


(A10) 


For the case of a constant surface temperature , both the thickness of 
the thermal and flow boundary layers are zero at the leading edge of the 
plate; thus, C equals zero in equation (AlO) and 


A 

8 



(All) 
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APPENDIX B 

VERIFICATION OF BOUNDARY CONDITION 


Equation (47) in the text is 


t(x,y) - By 1 / 7 


*° + 1 


, . . 39/280 —1/40 

t(L,0) x L dL 




4/35 10 J ( ^312/280 ^ ^273/280 ^39/280 ^ 43/39 


; 


where 


280/273 


U _ ^273/280 _ b 39/4 y 39/28 ^-39/280 ^ 

To simplify the integration let 

_ ^312/280 _ l 273/280 x 39/280^ y n 


/ 


(Bl) 


(B2) 


(B3) 


or 


S* = _ 213 yl / 9 / 2e ° ir 1 / 40 ^ 
280 


(B4) 


When the terms of equation (B4) are transposed 


z 39/280 L -l/40 dL 



dq> 


(B5) 


From equation (B3) it can he found that 


L = 


f 39/40 


-39/280 -n 

X ' y 


0 


40/39 


(B 6) 


With this transformation of the variable, equation (Bl) becomes 



28 


t(x,y) = 


By 1/7 t Q 

4/35 


39 


_ 39/4 n+39/28 

B y 


39/40 —39/280 ] 

[ — x y 9y 


40/39 


x 312/2SO y n 


43/39 — 43E/39 

9 y 


When n = —39/28, equation (B7) reduces to 


t(x,y) = 


By 1 -' 7 t, 

4/35 


39/4 


iffl 

39 


B 


312/280 

39/28 


v 40/39 

39/40 —39/280 39/28 ) 

x -x y <P/ ,<■ 


9 


43/39 


At y = 0, equation (B8) Becomes 


t(x,o) = r 

39 i 


39/4 


t(x,0)dcp 


9 


43/39 


On integration, equation (B9) becomes 


t (x,o) = t(x,o) j 22. cp^ /39 
39 \ 4 


39/4 

B 


which reduces to the identity 


t(x,0) = t(x,0) 
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0 y~ n+1 / 7 d(p 

(B7) 


(b8) 


(B9) 


(BIO) 


(Bll) 
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Figure See so* s original data , reference 
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Figure 2.- Scesa's data altered by approximate theory with m s 0. 
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Figure 3. - Scesa’s data altered by approximate theory with m 
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Figure 4.- Values of Pflt) defined by equation (34). 
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n 

Figure 6. - Numerical values of Y n defined b y 

equation (65). 
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Figure 7. - Effect of variable surface temperature on 
the heat - transfer characteristics in example 

problem. 
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